Charge and spin density response functions of the clean two-dimensional electron gas 
with Rashba spin-orbit coupling at finite momenta and frequencies 
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We analytically evaluate charge and spin density response functions of the clean two-dimensional 
electron gas with Rashba spin-orbit coupling at finite momenta and frequencies. On the basis 
of our exact expressions we discuss the accuracy of the long-wavelength and the quasiclassical 
approximations. We also derive the static limit of spin susceptibilities and demonstrate, in particular, 
how the Kohn-like anomalies in their derivatives are related to the spin-orbit modification of the 
Ruderman-Kittel-Kasuya-Yosida interaction. Taking into account screening and exchange effects of 
the Coulomb interaction, we describe the collective charge and spin density excitation modes which 
appear to be coupled due to nonvanishing spin-charge response function. 
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I. INTRODUCTION 

One of the working principles of semiconductor 
spintronica«»a is based on the idea to exploit spin-orbit 
(SO) coupling for a manipulation of an electron's spin by 
means of electric fields. The SO coupling of the Rashba 
type^ arises in a two-dimensional electron gas (2DEG) 
at semiconductor hetero junction due to the quantum 
well asymmetry in the perpendicular direction, and the 
strength an of this coupling can be tuned by a gate 
voltage^. 

A theoretical description of SO-related phenomena 
in the 2DEG is provided by coupled transport equa- 
tions for charge and spin components of the distribution 
function£££&i£. In the regime of a linear response to 
external fields these equations appear to be intimately 
linked to the density response functions such as charge 
and spin susceptibilities and - more peculiar - spin- 
charge response functions. In the presence of impurity 
scattering all of these functions have been previously 
evaluated in the quasiclassical approximation in the both 
diffusive^ (q < ^) and non-diffusive±i < q < 

&f) regimes, where q is a momentum transfer, Hf and 
vf are Fermi momentum and Fermi velocity, and r is an 
clastic scattering time. 

Recently it has been also remarkedii»i2> that the qua- 
siclassical results for the Rashba system are validated 
only in the presence of a finite amount of disorder such 
that t _1 ^> m*a 2 R , where m* is an effective electron's 
mass. Therefore, they cannot be straightforwardly ap- 
plied in the extreme collisionless limit r — > oo even at 
small q -C fcjr, and the response functions of the clean 
2DEG with Rashba SO coupling require a more refine 
consideration at finite values of q and frequency ui. 

Various response functions in the clean case r — > oo 
are most easily evaluated in the long-wavelength limit 
q -» q^u^jm On the other hand, the knowledge of 
the dynamic response functions at finite q enables one to 
find dispersions of the collective charge and spin density 
excitation a 17 ' 18 i 19 ' 20 i 21 ' 22 i 23 occurring in the presence of 



electron-electron interaction. In particular, in Ref. |23| 
the polarization operator, or the charge susceptibility, of 
the system in question has been calculated analytically at 
arbitrary momenta and frequencies, and the SO-induced 
attenuation of the charge density mode (plasmon) has 
been quantitatively described within the random phase 
approximation (RPA). 

In this paper we analytically evaluate (Sec. IIII[) the 
other density response functions following the compu- 
tational scheme elaborated in Ref. [23|. In Sec. IIVI we 
derive the static limit of the spin susceptibilities, and ob- 
serve an occurrence of the Kohn-like anomalies^ in their 
derivatives. We demonstrate how they are related to 
the SO modification 2 ^ of the Ruderman-Kittel-Kasuya- 
Yosida (RKKY) interactio n 26 ^ 7 ' 28 between two magnetic 
impurities. In Sec. [V] we compare our expressions for 
the dynamic response functions with the results of the 
long-wavelength and the quasiclassical approximations, 
and make conclusions about applicability ranges of the 
latter. Finally, in Sec. I VII we revisit the problem of 
the collective charge and spin density excitations treat- 
ing electron-electron interaction in terms of the Hubbard 
approximation 2 ^, i.e. our consideration extends beyond 
the RPA scheme. We demonstrate that the charge and 
the spin components are coupled in the obtained collec- 
tive modes, which is a consequence of a simultaneous ac- 
count of the non-zero spin-charge response function and 
the exchange vertex corrections. 



II. BASIC DEFINITIONS 

Let us consider the 2DEG with SO coupling of the 
Rashba type- which is described by the single-particle 
Hamiltonian 



H k = 



2m* 



auk h 



k • 



(1) 



where /ij. = a 



sin (/>k — a v cos 4>\s. is the spin-angular part 
of the Rashba spin-orbit coupling term, and we use the 
units such that K = 1. 
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The spectrum of (QJ is split into two subbands, 
A- 2 



k 2m* 

the corresponding eigenstates being 



1 



V>k- = —7= 



(2) 



(3) 



The matrix Z-/k diagonalizing the initial Hamiltonian (|TJ) 
as well as /i^ , i.e. a z — U^h^U^, is then given by 

Uk = ("0k+ ; vfac.- ) = --j= [1 - ^ cos cf>u - ia y sin k ] ■ (4) 

It is also convenient to introduce the projectors onto the 
eigenstates 



Pk± = tpk± ® V'ki = 



l±hl 



(5) 



which allow us, for example, to spectrally decompose the 
Hamiltonian = J^u=± e k^k^ as well as the (retarded) 
Green's function 



Gret 
k£ 



iO 



(6) 



A linear response of the charge density (a = 0) and 
the spin densities (a = x, y, z) to an external spatially in- 
homogeneous and nonstationary perturbation V^ u , which 
consists of a scalar potential ((3 = 0) and a magnetic field 
(/? = x,y,z), is usually determined in the framework of 
the Kubo formalism 29 . Applying the standard technique 
of the linear response theory and using for convenience 
the representation ©, one can establish an expression 
for the (retarded) density-density response functions, 



n,n'=± 



d 2 k n F {e£) - n F (e£ ) 



(27T) 2 ^ + , + ^ - £ ^ +q 



V 





k,k+q;/i,/x' ' 

(7) 



where 



(8) 



are the overlap functions; n F denotes the Fermi distribu- 
tion, and cr° = 1. 

In the explicit form , are listed in Appendix 

B of Ref. I 111 . Here we quote their symmetry property 
which can be directly established from the definition ||5J): 



■a/3 

k,k+q;/^,^A 



a/3 

— k— q, — k;/^',^ 1 



(9) 



where s = 1 for the charge-charge and spin-spin compo- 
nents, and s = — 1 for the spin-charge components. 

The expression ([7]) includes definitions of a polarization 
operator (a, ft = 0), spin susceptibilities (a,/3 = x,y,z) 
as well as of spin-charge response functions (a = and 



(3 = x,y,z, or vice versa). In the presence of SO cou- 
pling the latter functions do not vanish, and their study 
represents an especial interest. 

We note that the expression ([7|) can be alternatively 
found in terms of the equations of motion for the local 
charge and spin densities (see Appendix [A] for details). 
A matrix formulation of this approach provides a conve- 
nient tool for an account of screening and exchange ef- 
fects in presence of electron-electron interaction. In more 
detail this will be discussed in Sec. IVII 



III. EVALUATION OF 

The functions ([7]) have been previously treated at finite 
q in terms of different approximations. The most typical 
of them are: 1) the small-q (long-wavelength) formal ex- 
pansion of the whole integrand (see, e.g., Refs.l6lll3rk and 
2) the quasiclassical approximation (see, e.g., Ref. fl7l ) 
which is usually performed in the quantum kinetic equa- 
tion approach. 

In Ref. the polarization operator x°° in the clean 
limit has been evaluated beyond these approximations, 
and the obtained result has been used for an estimation 
of the accuracy of the long-wavelength expansion. In 
particular, the latter has been shown to be applicable in 
the limited range of the very small q -C k R /kp, where 
kp = y / 2m*ep + k 2 R and kp> — m* ccr is the Rashba mo- 
mentum splitting. 

As for the quasiclassical approximation, it has been ar- 
gued in Ref. YLVl that its application in the presence of SO 
coupling is validated at the finite values of the disorder 
broadening r" 1 3> k 2 R /m* , which smoothens the diver- 
gences of the quasiclassical result near the boundaries of 
the SO-induced (intersubband) particle-hole excitation 
region in the (q, aj)-plane. One still might hope that the 
quasiclassical approximation is trustful in the extreme 
collisionless limit r — > oo, provided one does not come 
too close to the boundaries in question. For this reason 
we are going to revisit its accuracy in the context of our 
present calculations. 

It also remains unclear how the two above mentioned 
approximations are related to each other in the clean 
limit. Both of them are elaborated for small values of q, 
but seem to give different results even at q — > 0. For ex- 
ample, an application of the long-wavelength expansion 
to the optical conductivity yields a box-like functio n 13 i 14 
which is finite at frequencies 2afjkp — 2m*a 2 R < u) < 
2a^kp + 2m*a R , while in the quasiclassical approxima- 
tion the width of this frequency window cannot be re- 
solved at all. 

In this section we are going to evaluate the functions 
(UJ without making any kind of approximations. We 
will neatly follow the computational scheme elaborated 
in Ref. [23; for x°°- Thus, we will derive analytic expres- 
sions for the other response functions and discuss their 
limiting behavior in the subsequent sections. 
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We start from the observation that, due to the mo- 
mentum space isotropy of the spectrum @ , the response 
functions ([7|) can be represented in the form 



Y Q/3 = V Y a/3l/i 



rf 2 k , ... . i-« 

a/3 



(10) 

(11) 



k,k+q;/i,/i 



■a/3 

k,k+q;/i,— /i 



"k+q 



A(w + iO) 



m _ -a* 

f k £ k+q 



+ X(w + iO) 



where the factor A~2~ originates from the property (9j. 
It is convenient to choose the basis in the momentum 
space such that q is aligned with x-direction. Then the 
matrix \ becomes sparse, the non-vanishing terms being 
X 00 , X xx , X vv , X zz , X° v = X y °, and X zx = ~X XZ ■ Intro- 
ducing x = cos(<?!>k — q ) = cos (j>, we obtain the following 
expressions at zero temperature 



oo, n 

Xnuj ,\ 

ZZ,[1 

Xquj, A 



1 

8^ 



-. _i_ k+qx 
lzt Tk+qT 



kdk 



-i k+qx 
iT |k+qT 



(12) 



"k+q 



A(w + iO) 



f k f k+q 



A(w + iO) 



Xqlj, X 
XX, p, 

Xqu,\ 



1 ± 



1 

8^ jo 

fc(2a: 2 -l)+g3: 

[k+qj 



fcdfc 



2tt 





it 



fc(2z 2 



(13) 



|k+q| 



f k+q 



+ X(u + iO) 



"k+q 



X(u + iO) 



l Xqu,\ 



8tt 2 

q+kx 



-flk R 



|k+q| 



XT 



271 



dcj> (14) 



|k+qT 



k+q 



X(u) + iO) 



C k+q 



X(uj + iO) 



Note that the components of x m the arbitrary basis can 
be recovered by an orthogonal rotation in the x — y plane 
(see Appendix B of Ref. [llh . 

After simple algebra we eliminate the odd powers of 
|k+q| in Eqs. (fL ^ -(fl4 |) . It means that the corresponding- 
integrands happen to be rational functions of k and cos <f>. 
Let us also note the identity 



X 



X" 



X 



x~ 



(15) 



which allows us to express, say, \ xx m terms of the other 
diagonal components. There remain, in fact, only five 
independent functions x = {x°°i X 0y , X VV i X zz i ix zx }> 
which can be conveniently labeled by the index j — 
1,...,5. Like in Ref. |23| . we also introduce the index 
i = 1,2,3,4 which denotes different combinations of 
{fi, X} = {-, +}, {+, +}, {-, -}, {+, -}, respectively. 



Defining the dimensionless units y = kj^/kp, z — 
q/2kF 7 v = k/k Fl and w = m*u/2k%, we cast (|T2 ]) -(|T3 |) 
into the form 



-ImX 

v 



(j) _ 



vg < f\v,z,w,y)dv, (16) 



1 . 



-Kzx?' = /, 



vf^(v,z,w,y)dv,(17) 



is the density of states in 2DEG 
(?) (i) 

per each spin component. The functions g\ and j\ are 



where v = v^d 

per each 
given by 

XC U) f 2» 
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U) 



dip sign(2?jz/c — [iyv + 2(z 2 — Xw)) 



X(x + Sl j) ) 5{x 2 +/9iX + 7i), 



U) _ 



C" 



0') I-21T 



, d(/.-5 — --5 

27r ./n x 2 + fta; + 7, 



where the coefficients 

2(z 2 - Xw) - fiy(v + fiy) 



Pi = 



7» 



(z 2 — Auj) 2 — f-iyv(z 2 — Xw) — z 2 y 2 



(18) 
(19) 

(20) 
(21) 



are the same for each j, and S(. . .) in (fTS)) denotes the 
Dirac delta function. The difference between the re- 
sponse functions appears only in the form of the 
coefficients and which are listed in Ap- 

pendix [B] for all j's. 

We note that the real part of x can be represented 
as a sum 



Re X 0') = x u) + Rex UM + Rex 



(22) 



where the term x^ = ~ v ^2if^ 1S nonzero only for 
j = 3 and j = 5 [see Eq. (|BlT|) ]. The terms Rex (i) ' 7 
and Rex^' J/ are obtained by integrating the functions 
Ei and « / i °' ) " ri , where = ) e(/3f-4 7i ) 

and/«' // = /«6(4 7l -/? 2 ). 

Performing angular integration and a subsequent 
change of the variable v — ► r(w) according to Eq. (34) of 
Ref. [23, we obtain expressions for Imx^ and Rex^' 1 in 
the form of x^ found previously, 



-0*) 



-Imx 0) = V a / c2t£ W+ (t) 
+G(l-4w)^2 dTC {j) -{T) (23) 

,, — 1- Jt+-(u) 



r T +-iv) 



i-t+-(v) 

+2Q(y 2 -Aw) / dT£ (j) "(r) 



4 



and 



-Rex (i),I = V / d T K u)+ {T) 



+0(1 -Aw) / dTTZ {j) -{T) (24) 

/.T_ + (0) 

+2Q(y 2 - Aw) c1tTZ (j) -(t), 



where 



Tl,2 


= ±w/z, 


T3A = —y±z 


; 


(25) 


T (T A(a;) 


1 

~~ 2 


—x - 


f o"v x 2 + AXw 


) 


(26) 


0)± (r) 


= C^(t) 


sign(r 2 + yr ± 


w), 


(27) 




= K^(t) 


sign(rz =p u;(r 


+ y)/z)- 


(28) 



In the representation (|23|) - (|24|) the actual integration 
limits are universal for all response function. The differ- 
ence appears only in the form of the integrands 



K u) (t) = Q (j) (t) 



Q(-P(t)) 



4 

p(r) = n( T -^)' 



k=l 



which are specified for each response function by 
1 



2z 



(t -t 3 )(t -t 4 ), 



Q (1) (r) 
q( 2 )( t) = --QW (T ), 

ZT 

Q (3) W = ^Q«(t), 
Q (i) (r) = A(T-n)(r-7a) 



)(5) 



(r) = 



2r 2 



Q W (r). 



Interestingly, for the function x 
would obtain the term 



(29) 
(30) 

(31) 

(32) 
(33) 

(34) 
(35) 
(36) 

X'~ z one 



QW(r)~Q^(r)+Q (i Hr) = ^^Q^(r), 



(37) 



which is anticipated after comparison of Eqs. (I35|) and 
(HI with flMI- (EH- 

Let us make several comments about the obtained re- 
sults ([23]) and ([24]) . 

1) First of all, note that the overall sign in the second 
line of (|23|) differs from its counterpart in the correspond- 
ing equation (35) of Ref. [23j. We use the present oppor- 
tunity to correct the misprint in the previously derived 



expression. Fortunately, it did not affect any other result 
of Ref. [H. 

2) The actual intervals of integration in Eq. (|23|) are 
explicitly written down in Eqs. (1B23P - (|B26[) . None of 
them contains the point r = 0, which means that one 
should not worry about the convergence of integrals over 
£(3,4,5) ^ ^ x/r 2 near this point. 

3) The explicit analytic relations for Im%W in terms 
of elliptic functions^ can be found for all j's in the same 
fashion as it has been done before for x [see Appendix 
C of Ref. S3. 

4) Some of the actual integration intervals in Eq. (pM]) 
do contain the point r = 0, which means that the cor- 
responding integrals J drTZ^' 4 ' 5 ^ (t) are divergent in its 
vicinity. However, the whole expression Eq. (|24|) is con- 
vergent and well-defined, since the singularities exactly 
cancel each other. In order to make Eq. practically 
useful, one has to substitute J drlZ^ (t) by the corre- 
sponding difference of primitives, which can be also found 
in terms of elliptic functions. 

5) Eq. contains only the contribution Kex^' 1 
to the full function Re%w. The contribution x is 
quoted in (|B11[) . and it remains to calculate the contri- 
bution Rex^' 11 ■ Making a complex change of variables 
t = + fiy) + iy/Aw- (i> + ny) 2 } (cf. Eq. (34) 
of Ref. 23) one can as well find Kex^' 11 in the analytic 
form which would involve the same integrand TZ^' (r) and 
a path of integration lying in the complex r-plane. Omit- 
ting technical details of this evaluation, we present the 
explicit expressions for Rex^'' 7 +Kex^'' 11 in Eqs. (|B29j) - 
(|Bl2|) . 



IV. STATIC LIMIT 

In order to find a static limit of x Q/3 one should consider 
with caution Eqs. (|B29j) and (|B32I) at w -> 0. One can 
then find 

- v~ x lim x 00 = 2 + -0(2/ - \z- l|)sin^ 

w— *o 2 

— ^] &(z — (1 — ny)) (/i?/v smi/j + cos^ + 21^ cos ip) 
n=± 

— 2Q(z — 1) cos^ arccoshz, (38) 
v- 1 lim X zz = 2 + ^ 9 ( z ~ (1 - 



cos ip 



At=± 



+ 2©(2__L) / arccos i lz _ ^1 _ (l/z) 2 ") 
cosw V " / 



cos-0 



(39) 



v- 1 lim x yy = 2 - 20(z - 1) cosi/Vl - (l/^) 2 , (40) 



ii/ 1 lim y za: = ^r0(y- \z-l\ 

w — *0 Z 



Y @{z- (1- W)) (m^V - 2Z M tan-0) 



+20(z - 1) tan V (arccoshz - ^/l - (l/z) 2 ) , (41) 
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FIG. 1: Diagonal components of the static spin susceptibility 
near z = q/2kF ~ 1 plotted at y = fefj/fci? = 0.1. Solid, 
dashed-dotted, and dashed lines correspond to a = x, a = y, 
and a = z, respectively. The inset shows the off-diagonal 
component — i\ zx {z, 0)/v. 



where sixiip = y/z (for y < z), sin-0^ = (1 — [iy)jz 
(for 1 — fxy < z; note that this notation differs from its 
counterpart in Ref. 23= by /i — > —fx), and 



^ = l^(z) = In 



1 + zsin(V> M + ^-0) 
2V2z cos | cos ^0 



(42) 



The static limit of can be found from the identity 
(ITS]) . The off-diagonal spin-charge term x 0i/ identically 
vanishes in this limit, which means a decoupling of charge 
and spin components at zero frequency. 

On the basis of the derived expressions (f38|> - (|4T|) one 
can observe that for z < 1 — y all diagonal terms are equal 
to x aa — — 2^, while their large-z asymptotes are x™ ~ 
—v i . The off-diagonal spin-spin term \ zx equals zero 

at z < 1 - y, and \ zx ~ ^ (l + ^p) at large z » 
1. The behavior of all components of the static spin 
susceptibility near z ~ 1 is shown in Fig. [TJ One can see 
that x vv has a discontinuous derivative at z = 1, while 
the derivatives of x xx \ X Z1 ' ■> and X zx are discontinuous at 
z = 1 =F y. These anomalies are analogous to the Kohn 
anomaly of the polarization operato r 23 ' 24 ! 29 at z = 1. 

Using ([39 ]) -(|4T |) we can find a SO-modification 
of the Ruderman-Kittel-Kasuya-Yosida (RKKY) 
Hamiltonia n 26 i 27 i 28 , which describes an indirect ex- 
change interaction between two localized magnetic 
impurities. In general case the RKKY Hamiltonian 
reads 25 



ttRKKY 
#1,2 



jrkky 

a./3—x,i 



(43) 



where Si j2 are the spin operators of impurities, and i"i2 = 
ri — r 2 is the distance between them. 

In the presence of SO coupling the Hamiltonian (14"3")) 
becomes anisotropic in spin space, since the matrix \ 



is no longer proportional to the unit matrix. Let us 
find asymptotic values of X a ^( r i2) at large ri 2 3> kp 1 . 
For simplicity we assume that the vector ri2 is aligned 
with x-direction in the coordinate space. Inspecting (|39p ~ 
(|4Tj) and restoring the dimensional units, we establish the 
asymptotic form of the right-sided derivatives 



d X yy 






1 2(? c0 




dq 




3 2k F \ 


q - Qco ' 




dx zz 




d\ xx 




/ 2g CM 


dq 




dq 




/ q - q cft ' 


dx zx 




ifiis 


/ 2g CAI 




dq 




3 


q - 9cp 





(44) 
(45) 
(46) 



near the discontinuity points q C Q = 2k F and q cfi — 2k F — 
2/ikfi. Performing the Fourier transformation of x°^(q)t 
we obtain the following leading asymptotic terms 



X VV ( r 12) 

X zz (ri2) 
X ZX (ri2) 



v sin(2fci?ri2) 



(47) 



xx( . i> V- gc M sin(g c;j r 12 ) 

£t=± ^ 



sin(2fcF^i2) 



cos(2fc fl ri 2 ), 



9c M cos(g CAt ri 2 ) 

"2tt ^1 ^ 



2k F r\ 2 



v sin(2fcpri 2 ) 



sin(2fc fl ri 2 ), 



(48) 



(49) 



12 



which oscillate in the coordinate space with the periods 
27r/g c o and 2ir/q cfl . Substituting them into Eq. (|4"3")l . we 
obtain 



H 



RKKY 
1,2 



F 2 (r 12 ) 



E 

a,0—x,y,z 



s?o a P(e 12 )s%, (50) 



where the range function 



rp f \ 7 2 vsm(2k F r 12 ) 

F 2 (n 2 ) = -J RKKY 2 ( 51 ) 

~ r 12 



is the same as in the absence of SO coupling. The SO- 
modification of the Hamiltonian (I43|) consists in the spin 
twist determined by the orthogonal transformation 



cos 9 12 — sin 12 
()[f) l2 ) =[01 

sin 6*12 cos6>i 2 



(52) 



with the rotation angle # 12 = 2/cRr 12 . The expression 
(|50p is in agreement with the corresponding result of 
Ref. HI 
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V. BEHAVIOR OF THE RESPONSE 
FUNCTIONS x al3 AT SMALL MOMENTA 

A. Exact expressions for lmx al3 in the SO-induced 
particle-hole excitation region 

It has been discussed in Ref. that an account of the 
Rashba SO coupling leads to an extension of the bound- 
aries of a particle-hole continuum, or Landau damping 
region, which is defined by Im;^ 1 ) ^ 0. It has been also 
established that this extension has a shape of the wedge 
bounded by the parabolas — (z — y) 2 — (z — y) = Wi{z) < 
w < wi{z) = (z + y) 2 + (z + y) [see Fig. [2]. 

Since the representation (|12p - (fT4|) manifests the same 
pole structure for all response functions, their imaginary 
parts appear to be nonzero in the same domain where 
Im X ( 1 ) ^ 0. Analyzing ((B25|) . we extract an explicit 
expression for Imx^ in the SO-induced particle-hole ex- 
citation region 



1 

0.8 
0.6 
0.4 
0.2 



/'i 

J 
■r % 



,1 \ :\ 

< \ w 
v 



■■\\v s 



0.05 



0.1 



0.15 w 



FIG. 3: Spin-galvanic function (|54|l at y = 0.1 and the values 
of z = 0.0, 0.025, 0.005, 0.01, 0.02, 0.03, 0.04. In particu- 
lar, the solid line depicts the long- wavelength limit, and the 



dashed curve corresponds to z = y 



0.01 



1 



j m ^(j) — — 0(w 2 — w)Q(w 



V) 



-0(u;i - w)Q(w ~ w 2 )A l i} (-h) 

where 101,2 = ( z iy) 2 i (z±y) and 103,4 = — (z±y) 2 ±(z± 

y). The functions and the arguments t\ and t% are 
defined in Eqs. (|BT7|) . |B27)l and ([E28|) . respectively. On 
the analogy of A^} explicitly quoted in Ref. HH, one can 
as well express the rest A^} in terms of elliptic functions. 



(53) 



0.15 




0.05 



z* 0.02 



0.06 



0.1 z 



FIG. 2: SO-induced extension of the particle-hole excitation 
region at y — 0.1 - the light-gray area bounded by the parabo- 
las wi = w+ and W4 = W- . The small darkened triangle inside 
it indicates the range of applicability of the long-wavelength 
approximation. The dashed lines depict quasiclassical bound- 
aries w± = y±z. 



For example, in Eqs. (|B18[) and (|B19|) we write down 



explicit formulas for A 2 ' and A)?' respectively. 



(4) 



Eq. (j53[) allows us to study the behavior of the 
density-density response functions at finite momenta. On 
its basis we can also describe the spin-density response 



2M}j.E x to a longitudinal electric field E x 



dx 



(or E x = —iqV° in the momentum representation). The 
spin-galvanic response functio n 31 i 32 i 33 M% is then related 
to x° v via M x = iex° v I '(2<z) ■ Let us also introduce the 
rescaled function 



16A* 



ReM^ 



Mi 



Im% 



0y 



(54) 



Using (|B18|) . we calculate and plot it in Fig. [3] at fixed 
finite values of z. We remark that its frequency depen- 
dence in the range 104 < w < w\ corresponding to in- 
tersubband transitions is similar to that of the fmite-q 
conductivity studied in Ref. [231 . 



B. Long- wavelength limit 

Spatially uniform spin susceptibilities of the 2DEG 
with SO coupling have been previously considered in 
Ref. [TH We recover the corresponding expressions in 
the long-wavelength limit 



X 

lim = —2 lim 

z^0 v z^0 



X 



•x(yy) 



where 

r(w) = In 



{w-y 2 ) 2 



{w + y 2 ) 2 -y 



= 2 +2^M, (55) 



+ iirQ(y 2 - \w-y\). (56) 



The small-g behavior of the polarization operator of 
2DEG with Rashba SO coupling has been approximated 
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in Refs. HIH by the expression 



— ~ — 9 + :r~ u^h 



(57) 



which is obtained after the formal expansion of x°° m a 
series of z = <C 1. Later on it has been remarked 2 ^ 
that the formula (|57|) is, in fact, reliable only for the 
values z < z* « y 2 , where z* is the point of intersection of 
the parabolas 102(2) and u>3(z) [see Fig. [5]. Analogously, 
we can find approximate relations for the off-diagonal 
terms 



and 



zw 



Ay 



r(w) + 



wy{n - 2y) 



z 

4y 



r(w) H 



(58) 



(59) 



which are also applicable in the small triangular region 
located at z < z* and bounded by W3 < w < 11)2- In 
particular, Eq. ([59]) accounts for the box-like shape of 
the function m% (JSJ) in the linit z — > which can be 
seen in Fig. [3l 



C. Quasiclassical approximation 



On the basis of the exact result (|53j) we can also esti- 
mate how accurate the quasiclassical approximation ap- 
pears to be, when it is applied to a description of the 
clean 2DEG with Rashba SO coupling. 

The quasiclassical approximation relies on the fact that 
all energy scales in the system are much smaller than the 
Fermi energy: qkp/m* ,u>,aRkp <§; ep. This inequality 
enables one to treat (| 1 2 1) - (fl^j) linearizing the branches 
of the spectrum near the corresponding Fermi points 
hp = kp — i±kpi and expanding ^ + , in a series of 
q. We note that a more systematic procedure of making 
the quasiclassical approximation is based on the gradient 
expansion of the quantum kinetic equation&£>^, which is 
applicable in more general - nonequilibrium - situations. 

Quasiclassical response functions x"' 3 01 the disor- 
dered 2DEG with Rashba SO coupling have been ex- 
plicitly calculated in Ref. [ll|. It has been also argued 
therein that the quasiclassical approximation is justified 
in the presence of a finite amount of impurities such that 
t _1 > k 2 R /m* . This condition confines r _1 from below. 
Had we ignored it, we would have obtained in the extreme 
collisionless limit r — > 00 the quasiclassical expression 
for, say, out-of-plane component of the spin susceptibil- 
ity in the form 



-Imx, 



qc 



9(z 2 -{w- w) 2 ) 
yjz 2 -(w- [iy) 2 



(60) 




0.04 0.06 0.08 0.1 0.12 0.14 0.16 w 



FIG. 4: Comparison of the exact (solid line) and the quasi- 
classical (dashed line) functions — Imx 22 /^ at z = 0.04 and 
y — 0.1. Vertical dotted lines correspond (from the left to 
the right) to the frequency values W4, 11)2, 11)3, and w\. Ver- 



tical dashed lines correspond to the values w_ 



and 



toi c = y + z demarcating the quasiclassical boundaries of the 
SO-induced particle-hole excitation region. 



The lines w± represent the quasiclassical boundaries 
of the SO-induced (intersubband) particle- hole excitation 
region, and they differ from the actual parabolic ones 
w + = w\ and W- = W4. In Fig. [2] w± are depicted by 



the dashed lines, and one can observe that 



lies 



which is divergent near w 



\v±z\ 



between wi and W3, and lies in between W2 and W4. 
It is also evident that in the quasiclassical approximation 
the finite basis y — y 2 < w < y + y 2 of the wedge at 
z = is not resolved. This means that the imaginary part 
of a quasiclassical counterpart of r(w) (|56p appears to 
be delta-peaked, and therefore the corresponding finite- 
valued results of Sec. IV Bl can not be reproduced in the 
quasiclassical approximation. 

It can be anticipated that the quasiclassical approxi- 
mation is still reliable in the triangular area u>2 < w < W3 
at z > z* « y 2 . In Fig. 0] we confirm this surmise compar- 
ing the function (|60| with the exact expression calculated 
on the basis of Eqs. ([53} and (|B19j) . 

Thus, we conclude that the long- wavelength and the 
quasiclassical approximations do not have any correspon- 
dence between each other, since they are applicable in the 
domains which do not overlap, i.e. at z < z* and z > z*, 
respectively. 



VI. COLLECTIVE CHARGE AND SPIN 
DENSITY EXCITATIONS 

Let us now consider the renormalization of the ma- 
trix x due to electron-electron interaction. Treating the 
latter in the Hubbard approximation^, we take into ac- 
count screening and exchange effects. An effective re- 
sponse matrix x [see Eq. (|A26|) ] is derived in Appendix 
lAl using the method of the equations of motion^. 
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In the basis specified by the condition q| le^ the matrix 
X can be decomposed into two 2 x 2-blocks 



X{Oy) 



x 0Q x 0y \ v _ (x xx x xz , „.,. 

-^yO yVV J ' A(a:z) I ^zx ^zz j ■ K^l) 



It follows from (|A26|) that the spin-charge X(0y) ancl the 
spin-spin Xixz) blocks are renormalized independently of 
each other. In particular, 



X(0y) 



1 



X(0y)^(0y)j X(Oy), 

X(xz) = y- + JX(xz)) X(xz), 



(62) 
(63) 



where F( 0y ) = diag{w 9 - J, -J} w diag{v 9 , - J}, v q = 

n 2 2 

is the Coulomb interaction, and J = -P- is the 
Hubbard's vertex exchange term at small q <C &f. 

Let us first consider the spin-charge block and find the 
renormalized response functions 



X 00 + JA {0y) 



Hov) 



X vv ~ v q A (0y) 



Hoy) 



~0y _ _X 



Hoy) 



where 



A 



(Oy) 



det X{ oy)=X 00 X yy -(x° y y 



(64) 
(65) 
(66) 

(67) 



A (0y) = det (1 - X(oy)F( y)) 



w (1 - v qX 00 ) (1 + JX VV ) + Jv q (X 0v ) 2 ■ (68) 
For the spin-spin block we have 

~x* _ X XX + J\zx) 



A 



X = 



X zz + JA 



(zx) 



A 



(zx) 



Hxz) 



where 



A ( xz) = detx ( «)=X x V + (x""r 



(69) 
(70) 
(71) 

(72) 



A 



(xz) 



det (1 + Jx(xz 



= (1 + J X XX )(1 + JX ZZ ) + J 2 (X XZ ) 2 - (73) 



Dispersions of the collective charge and spin density 
excitations are determined from the equations A(oy) = 
and A( xz ) = 0, and the response functions x a/3 are 
strongly enhanced at the parameter values satisfying 
these conditions. 

Before quantifying x Q/3 , let us qualitatively discuss the 
role of the spin-charge mixing term x° v as well as the role 
of the exchange corrections. Note that if either x° v — or 




.04 



0.02 



0.02 0.04 0.06 0.06 0.1 0.12 0.14 



FIG. 5: Contour plot of —\mx xx (z,w)/v. Parameters: y 
0.1 and r s = 0.6. 



J = 0, the dynamics of the charge and the spin-y densi- 
ties becomes decoupled, and the corresponding plasmon 
and SDEy (spin-?/ density excitation) modes are inde- 
pendently found from the conditions 1 — v q x 00 = and 
1 + Jx yy — (the latter equation makes sense at J ^ 
only). The role of x° y , whatever small it might be, is con- 
siderably strengthened when the dispersions of the plas- 
mon and the SDE y modes come close to a degeneracy 
point. In fact, x° v lifts this degeneracy, thus mak- 
ing a possible crossing of these modes avoidable. There- 
fore, we expect that an account of the spin-charge mixing 
along with the exchange interaction would lead to non- 
trivial features in profiles of x 00 and x yy near the avoided 
intersection. It is also implied that the peaks correspond- 
ing to the both collective modes would appear in every 
response function of the spin-charge block. In particular, 
the spin susceptibility x vv is expected to manifest reso- 
nant features at the position of the plasmon dispersion, 
while x 00 should have a peak corresponding to the SDE y 
mode. 

We remark that the decoupled plasmon mode in 2DEG 
with Rashba SO coupling has been previously considered 
in Refs. [19.20.21.221231 in the random phase approxima- 
tion (RPA), i.e. at J = 0. In turn, the coupled spin- 
as - spin-z collective modes, which are determined by 
A(xz) — 0; occur at J ^ 0, i.e. their description re- 
quires an extension of the RPA. Their dispersions, as 
well as the dispersion of the decoupled spin-jy mode (at 
neglected spin-charge mixing x° y — 0), have been pre- 
viously considered in Ref. U3 in terms of the Hubbard's 
approximation with the bare spin susceptibilities calcu- 
lated in the quasiclassical approximation (see Sec. IV C) . 

Using our exact expressions for x"^; we are able to 
study dispersions of the collective charge and spin den- 
sity modes in more detail. They can be visualized, for 
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0.02 0.04 0.06 O.OS 0.1 0.12 0.14 Z 



FIG. 6: Contour plot of — lmx vy (z, w)/v. Parameters: y = 
0.1 and r s = 0.6. 



example, in contour plots of the response functions x a ^ ■ 
Absolute values of x a ^ providing a useful information for 
inelastic Raman scattering 34 are available as well. 

Let us start our consideration from the renormalizcd 
response functions (|69p ~ (|71[) constituting the spin-spin 
block X(xz) ■ In Fig- El we present the contour plot of 
— Imx xx /v at y = 0.1 and the Wigner-Seitz parameter 
r s = ^" e2 = 0.6 (note that J = ^^). One can ob- 
serve the dispersions of the two coupled spin- a; - spin-z 
collective modes. We state that their spectra are in the 
qualitative agreement with those previously predicted in 
Ref. [I?] on the basis of the quasiclassical approximation. 
However, the absolute values of x a ^ differ from their qua- 
siclassical counterparts Xqc ■ The reason is the same as 
discussed in Sec. IV CI for the case of bare susceptiblities. 

In turn, the response functions (|64p ~ ([M)) of the spin- 
charge block X(0y) manifest novel qualitative features due 
to the account of the spin-charge mixing x° y along with 
the exchange interaction. In Fig. [6] we present the con- 
tour plot of the function — Imx yy /V at y = 0.1 and 
r s = 0.6 (note that v q = 2 ^ zv ) • It also contains 
the two collective modes: plasmon-like and SDE y -like. 
The plasmon-like mode has almost the same dispersion 

w ~ \J 2^ as ^ e pl asmon rnode in the absence of SO 
coupling. The SDEj,-like mode originates at z = from 
the finite frequency value slightly below the bottom of 
the wedge. Being undamped, these modes do not come 
close to each other. They collide soon after the plasmon- 
like mode enters into the SO-induced damping region, 
i.e. at z > 0.023. In the vicinity of this point the spin- 
charge mixing x° y acquires its importance. Although a 
pictorial description of the avoided crossing loses its ob- 
viousness because of the modes' broadening, we prove 




0.07 0.08 0.09 0.1 0.11 w 

FIG. 7: The cross section of -lmx vy (z, w)/v at z = 0.027 
(solid line). The dashed line corresponds to the case of inten- 
tionally neglected spin-charge mixing term x° v ■ Parameters: 
y = 0.1 and r s = 0.6. 

that it does happen. For this purpose we plot in Fig. [7J 
the cross-section of — Imp^/f at z = 0.027. Thereby we 
show that instead of the only SDE y peak (dashed line) 
occurring at intentionally neglected x 0y = 0, we obtain 
at x° v 7^ the two well-resolved peaks (solid line) corre- 
sponding to the plasmon-like and the SDE y -like modes. 
A formation of the dip between them is interpreted as 
an avoided crossing of the two broadened modes. A rele- 
vance of such interpretation becomes even more evident, 
if one would compare several subsequent cross sections of 
— Imx vy /v and — Imx 00 ///. 

As one can see from Eq. (I5tl)) . the spin-charge mixing 
term is renormalized by electron-electron interaction as 




0.08 0.09 0.1 0.11 0.12 w 

FIG. 8: The renormalized spin-galvanic function (|74p at z — 
0.01 and the parameters y — 0.1 and r s = 0.6 (solid line). The 
dashed line depicts its non-interacting counterpart (r s — 0) 
given by (JS3J). 
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well. Let us define the function 

ml = i^-Imx 0y , (74) 

TTZV 

which is an interacting counterpart of the bare spin- 
galvanic function m v x ([54")) . In Fig. [5] we compare both 
of them at z — 0.01 and y = 0.1 in the frequency window 
y — 2y 2 < w < y + 2y 2 . We observe that at these pa- 
rameters the major effect of renormalization consists in 
a considerable amplification of the spin-galvanic function 
([74)1 due to the exchange interaction. 

VII. SUMMARY 

We have derived zero-temperature analytic expressions 
for the charge and spin density response functions of the 
clean 2DEG with Rashba SO coupling at finite momenta 
and frequencies. We have studied their static and long- 
wavelength limits as well as established the applicability 
range of the quasiclassical approximation. In the static 
limit we have observed the Kohn-like anomalies in the 
spin susceptibilities and showed how they are related to 
the SO modification of the RKKY interaction. 

The renormalization of the response functions due to 
electron-electron interaction has been considered in the 
Hubbard's approximation. We have studied the collec- 
tive charge and spin density modes which appear to be 
coupled due to the nonvanishing spin-charge mixing term 
X° y ■ One of the important consequences of this coupling 
is the emergence of the plasmon-like peak in the spectrum 
of the renormalized spin susceptibility. 
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APPENDIX A: EQUATIONS OF MOTION 

Let us define the local charge (a = 0) and spin (a — 
x,y,z) density operators 

Pq = 2^kq> Pkq = C k+q , (Al) 

k 

and derive the equations of motion^ for the expectation 
values p£ q = (/3 kq ). In the Fourier representation they 
read 

^kq=<[/Skq,if]>. (A2) 

The full Hamiltonian H = Hq + -ff C xt consists of the parts 
describing a system 

Ho = J2 E ^',o (A3) 

k',/3 



and an external perturbation 

q',/9 

where V is the spatial volume. For the Rashba system 
([I]) we have — (e k , a^ky, —ank x , 0). Later on we will 
also add the interaction term H lrA to the Hamiltonian H . 
Using the identity 

[Pkq>Pk',o] = cjc^^Ck'^k+q.k' - 4,(r /3 Cr a C k+q 5 kk /(A5) 

we find 

([Pkq,^]) = E^ + q(4^^ k+q ) 



- ^^(cjyW+q). (A6) 



It is convenient to introduce the representation 

^kq-([Pkq,tf0]}=£</4 (A7) 



in terms of the matrix A with the elements 

-^kq = ^ + ek-6k+q, 

AxO aOx 

-^kq ~~ "^kq — VLRSLy, 

^kq = A^ = a R q Xl 

-4£ q = --^kq = ia R {2k x + q x ), 

A lq = -•^kq= ia fi( 2fc J/+%)> 

< = < = -4k q = < = 0. (A8) 

In order to find the commutator of pg q with H cxt ()A4I) . 
we use the identity 

[Pkq.P-q'] = 4^^c k+q _ q , - c k+q ,^a Q c k+q . (A9) 

From the whole sum over q' in (|A4[) we pick out the 
only term q' = q. This is known as the random phase 
approximation (RPA). Then, 



(A10) 

and after averaging we obtain 

([/3£ qJ ^ext])«^E B W' ( AU ) 



where 

< = <4^c k - c k+q ^a«c k+q ). (A12) 

In order to fulfil the averaging, we have to transform 
c k = W k 7 k into the diagonal basis 7 k ± such that 

(7ku7V) = V/k 1 - (A13) 
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The transformation matrix lA^ is denned in Eq. ((4]) . Since 
the external perturbation is assumed to be small, we aver- 
age in (|A13|) with respect to the system's density matrix. 
Therefore we can identify / k = n f(e k ). 

The components (|A12|) are then found to be 



< = -^/ k ^Tr[(l + ^)^ ( TV^ k 



~ E /k +q Tr [(1 + ^)Ul + ya^ +ti 



- \ E ^ +q Tr K 1 + . ( A14 ) 

or, more explicitly, 

Skq = = E ^ (/£ Sin <^ - /k\q M 



S kq = S kq = - E * (/k C ° S ^ k ~ -^k+q COS 'Ak+qj , 

^kq = -#k q = t E i" (/k cos k + / k \ q cos k+q ) , 



BP = -52 = 



^kq 



kq 



•Em(/, 



°kq - °kq - °kq ~ °kq ~ U 



Q sm (p k + JQ +(i sin </> k+q ! 

(A15) 



Combining (|A7]> and (jAllj) . we derive the following 
matrix equation 



AcqPkq — — BkqVq, 



(A16) 



where the upper indices are omitted for brevity. Inverting 
the matrix „4k q and summing over k, we obtain p q = 
(Pq) = xKj, where 



11 E^kq^kq 



V 



(A17) 



In the mean field approximation we obtain 

([Pkq.tfiat]) = y EVPS' X 

q' 

x {(c k cr Q c k+q _ q -> - (c k+q ,a a c k+q )} 

+ ^ E;o'{(4 +q '^ c p')(i- q ^^ tt ^) 

(c k a^%,)(4' + q^ c k +q+ q'}} ■ (A19) 



q',p',/3 



The first sum in (|A19|) corresponds to the direct 
Coulomb term. We treat it further in the RPA picking 
the only term q' = q out of the whole sum. We obtain 
the following contribution 



([Pkq'-ffir 



\RPA 



(A20) 



— |(c k cr c k ) - (c k+q cr Ck+ q )| = ^-O kq , 

which accounts for the effect of screening. 

From the second sum in (|A19[) we can extract the ex- 
change self-energy term and the Hubbard's exchange cor- 
rection to the RPA. 

A contribution associated with the self-energy is ob- 
tained from (|A19j) after picking out the summands with 
p' = k + q' and p' = k + q, i.e. 

(LSkq,ff in t]> S = 

= ^?EV {(Ck+q'^k+q'KcX^Ck+q) 



q',/3 



E 



(4^ Ck+q )(c k+q+q ,^ Ck+q+q ,)} (A21) 
{S^(4^^c k+q ) - S^ +q (c k a Q ^c k+q )} , 



where 



yP - 



E V-kPq'O 

q' 



^ E V-k/^Tr[(l + //^Vl- (A22) 



The self-energy S k modifies the single-particle Hamilto- 
nian — ► + S k . After diagonalization we obtain the 
renormalized eigenvalues e k — > e k 4- S k , 



is a density response matrix of the non-interacting sys- 
tem. After the straightforward calculation, we recover 
from (|A17p the expression for the components of x- 

Let us now take into account electron-electron interac- 
tion 



Hint = ^ E E^'+q'^k'-q'^k' 



(A18) 



p'k'q' as 



E k = E V-k/q' [1 + MM'cos(0 k - ^y)], (A23) 
q':P' 

which correspond to the same eigenstates (j3|). Besides the 
shift of a chemical potential, the spectrum renormaliza- 
tion results in an effective value of the Rashba splitting, 
which has been previously studied in Ref. l35l . In our con- 
sideration we will, however, neglect this effect and discard 
the contribution (|A21|1 . 
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The Hubbard's exchange term is given by the sum- 
mands in (|A19[) with p' = k + q + q' and p' = k, i.e. 

<[/3£ q ,#int]) HUb = 



2V E W 1' { (4+q'^ C k+q+q') (<i +q tT a a C k+ ^ 

q',/9 

\ " nap 

7^ 2^ y k'-k^P k , q 



k',/3 



k',/3 



2V 



where u q is approximately regarded at small q as a con- 
stant = 2 J = 2|id. 



Collecting the contributions ()A20[) and ()A24[) and 
adding them to (|A16[) . we obtain the equation 



AiqPkq = rj #kq (Vq + Fp q ) 



(A25) 



where i* 1 = diag{t> q — J, — J, — J, — J}. Solving it, we find 
that pq = xVq, where 



X= {l-xF)-\ 



(A26) 



is a density response matrix of the interacting system in 
the Hubbard's approximation. 



APPENDIX B: EXPLICIT EXPRESSIONS FOR 
THE FUNCTIONS X U) 



Let us list the functions f- * , and 5^ which deter- 
mine the response functions x^' m Eqs. (fT5)l and (fTi?)) : 



/ (1 



1(2 
J i 

*? 

1(3 
•t i 

<? 
if 

if 



0, C 



(i) _ « - Mj/ 
2<; 2 z ' 



(z 2 — Aw) — /ij/u 
(v - py)z 



(Bl) 



A(ju - y) ( 2 ) _ fi\(Xw + [3izv + fiyv - z 2 ) 
; ° 5 ; — ~~ 



2z 



2v 2 z 



2 r 2 



(7if + ny)z 



Aw + flizv + jiyv — z 2 



y(v - A 1 ) ^(3) _ - /iyz + Pmyv 



2z 2 

z 2 — Aw + ji^,yv 
vz — [lyz + Pifiyv ' 

z 2 — Aw 



2v 2 z 2 



(w + ^iy)z' 
-2/ M 

(thv - yy)z 

Aw + (3izv — z 2 



ji-y r (5) _ /i(Aw + pjZV - z 2 ) 
2z ' 1 2v 2 z 2 



(B2) 



(B3) 



(B4) 



(B5) 




FIG. 9: The domains A, B, C, and D corresponding to the 
different orderings (fB6) - (fB9)l of the roots r k (f25l) . 



Employing (lBljl - (lB5j) in the framework of the compu- 
tational scheme elaborated in Ref. [H, one can derive 
Eqs. P5 ]) and f23 j) . 

The actual limits of integration in the latter expres- 
sions require further detailing, since the corresponding 
integrands £^(r) (29]) and ftW(r) & are defined at 
P{t) > and P(r) < 0, respectively. The polynomial 
P(t) (|3~1]) has the roots (|25p. which can be ordered 
differently depending on the values of w, z, and y. 

We identify the domains A, B, C, and T> in the plane 
(z,w) [see Fig. |2] such that 



A 
B 
C 
V 



r 4 < r 2 < n < r 3 , 

r 4 < t 2 < t 3 < Tl , 

r 2 < t 4 < t 3 < Tl , 

r 4 < t 3 < t 2 < Tl . 



(B6) 
(B7) 
(B8) 
(B9) 



We also define the unit-step functions 6 (-4) = 9(z 2 — 
yz - w), Q(B) = 9(w - z 2 + yz)<d(z 2 + yz - w)9(w + 
z 2 - yz), 9(C) = 9(w - z 2 - yz), and 6(2?) = (yz - z 2 - 
w), which realize projections onto these domains. Using 
these definitions we make the following decomposition 

x Ci) = x « + 9(^l)xi j) + 9(S)xi j) 

+ e(C) X «>+e(2%jf>, (BIO) 

where x'" 7 '' is nonzero only for j '■ = 3 and j = 5: 

X {3) = -2vy 2 jz 2 , X (5) = 2^y/z. (Bll) 

Inside of each domain the root's ordering is fixed, and 
one has to find out how T a \(. . .) occurring in Eqs. (|2"3"]) 
and (|24[) are arranged among t\ , . . . , T4 . 

Let us introduce x 4 < x 3 < x 2 < x\, where Xk's are 
identified with r^'s differently in each domain according 
to (|B6[) - (|B9[) . and define the primitives for the imaginary 



.4 



(x) 



dx'C^ix' 



X\ < x. 



(B12) 



A j) (x) 



j dx'C u) {x'), x 3 < x < x 2 , (B13) 
dx'C (j) (x r ), x< Xi , (B14) 
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and the real parts 

B[ j \x) = / dx"Jl^(x'), x 2 < x <x x , (B15) 

J X2 

B%\x) = / dx'K U) (x'), x A <x<x 3 . (B16) 

They have to be further detailed in each domain as well. 
For example, in the domain C we have 2; 3 = r 4 = — z — y, 
and therefore 



Ag(x)^A%(x) 



dx'C^(x'). 



(B17) 



After such a specification the primitives (|B12j) - (|B16j) can 
be explicitly found in terms of elliptic functions [see Ap- 
pendix C of Ref. HH . In particular, we quote the 
expressions for j — 2 and j = 4: 



k c [(w/z - y) 2 - z 2 



4zwy/w 
-(w/z) 2 n 2c Tl((p2c(x),n2c, k c ) 



{2wF (cp 2c (x), k c ) 



4^) 



where 



2(z 2 — y 2 )x V zx + w 
k c ({w/z-y) 2 ~z 2 ) 



[{x + y) 2 — z 2 



k c (z 2 - y 2 ) 



{zF{Lp 2c {x), k c ) 
kc)} 

E{<p 2c (x),k c ), (B19) 



4^/w(z 2 — y 2 ) 
- yn2jl(<P2c(x),n2c,k c )} 
zJw 



, \ ■ / (x + y + z)(w/z — y + z) 
ip2c{x) = arcsinW^ n / -,( B2 °) 



2(zx + w) 



k r . — 



ri2c 



2^B 



yj(z + w/z) 2 ~ y 2 ' 
2z 



w 



w/z + z — y' 



n 2c 



z(z + y) 



(B21) 



n 2c .(B22) 



(z 2 - y 2 )n 2c n.(cp 2c (x),n 2c ,k c )} ,(B18) Imx (i) in every domain 
I 



Let us now establish the actual limits of integration in 
Eq. J23J. Rewriting it in terms of l[BT2]) - (|Bl4l) . we obtain 



-Im X « = -e(w/z - tO [a£(w/z) - Agfa) + A«(-t0] - 0(i a - z - v)A$-(-t 2 ) 
+e(t 2 - z + y)Ag + fa) - 9(1 - Aw) [e(t 4 - z + y)A^ + (t 4 ) - Bfa - z + y)A$ + fa) 
- S(w/z - t 3 ) [a£(w/z) - Agfa) + Ag(-t 3 )] - Qfa - z y)A^-(-U)} , 



(B23) 



-Im X « = -B(z -y-h) [A$> (z~y)~ A$> fa)\ - Q(z y + h)A$ (-h) 
-(->(/, - : - h)Ag\z — y) — Qfa - z - y)Ag ) -(-T 2 ) + A^+fa) - 6(1 - Aw) x 



x {-0(t 4 - z - y)^ )_ (-t4) + ®fa - iu/«Mft )+ (*4) - 9(t 3 - iy/z)4i )+ (*3) 
-6(z - y - T 3 ) [4£(z — y) — Agfa)] + 0(w/z - t 3 )G(z -y + t 3 ) [a { £{z - y) 
-Q(w/z - U) [A 2 \\z -y)- A$j(-uj\ - Qfa - w/z)e{w/z - t 3 )Ag(z - y)} 



(B24) 



-Imx 



U) - Q, 



Qfa- w/z) Ag + fa)-A%-(-t 2 ) -e(z-y-h) A% (z - y) - A% (h) 



1(3)- 



-e(z - j/ + toe(z + y - *o4c (-*0 - @ (y - z - ^) A 2^ z - ») - ©C 1 - x 

x {e(t 4 - w/z) [^^ )+ (i 4 ) - A^-(-U)} ~&{z-y- h) [Ag(z -y)- Agfa) 
+0(z + y~ t 3 )Q(z -y + t 3 ) [4c (2 - 2/) - 4c Ha)] - 6(z + y - t 4 )6(z - y + i 4 ) x 
x 4^(z -y)- Ag(-U)] - Qfa - z - y)Q(z + y - t 3 )A ( g(z - y)} , 



0), 



(B25) 



-Im X « = A^ + {t 2 ) - Ag-(-t 2 ) + Ag~{-U) - Ag + fa). 



(B26) 
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In the above expressions we have used the definitions 

ti,2 = \ [VTTAw~ T 1] , (B27) 
t 3A = ^ [1 T Vl - 4wJ . (B28) 



The actual limits of integration for the real part (|24|) 
can be found in a similar way. However, Eq. (|24p accounts 
only the term Rex^ 7 . Complementing it by Rex*- 7 - 1 ' 77 
[see the remark 5) in the end of Sec. IHIj . we present the 
sum of the both terms in every domain 



-Re X « = 0(i 2 - z + y)B$(z - y) + 9(z - y - t 2 )B$(t 2 ) + &{h - w/z)B$( tl ) - 1B^{k) 



-[2 + B(w/z - tl )]B^(-w/z) - 9(ti - w/z)B%{-T x ) + 2B)j>{-i 2 ) - 0(z + y - t 2 )B%(-t 2 ) 
+29(y 2 - Aw) [l?«(t 3 ) + B$>(hj\ + 29(4™ - y 2 ) + fl^)] + 9(1 - Aw) [q(z - y - U)B<£{U) 

+ G(t 3 - w/z)e(z -y- t 3 )B$ (t 3 ) + [9(i 3 - z + y) + 9(i 4 - z + y)}B$ ( z -y)-G(z + y- t A )B$ (-t 4 ) 
- 0(t 3 - w/z)B^(-t s ) - 9( W /z - * 3 )fltj) (-«;/«)} + 9(4u, - 1) {^(fg) + S^(^) - b£>(-*§) 

>HS)} 



0')/ 



iCj)/ 



B 



U), 

2 a 



K\[z r 



z^J z 2 — y 2 



2s 



0) 



B${z-y)-B%(-w/z) 



(i), 



(B29) 



-B$(w/z) + 2B$(h) - 9(<r - z + y)B ( ^{t 1 ) - 9(y - 



t!)B$ (-h) -e(z + y- t 2 )B%\-t 2 ) 



+2Bg\-t 2 ) - 0(w/z - ti)B$>(-w/z) - 0(ti - w/z)Bg\-t 1 ) - 29(y 2 - Aw)e(w - z 2 ) [b$ (t 3 ) + B^ft) 
+29(y 2 - 4 W )9(z 2 - «,) [b«(* 3 ) + B^(i 4 ) ~ 2B£\-w/z)] - 2Q(Aw - y 2 ) [b$(%) + B^{t%) 
+9(1 - Aw) { [0(t 3 - to/z) + 9(* 4 - w/z)]B^ (w/z) + Q{w/z - t 3 )Q(t 3 - z + y)B { ^ (t 3 ) + Q{w/z - t A )B^ (t 4 ) 
+[9(t 3 - w/z) +Q(U- w/z)]B$ (-w/z) - 9(z + y - U)&{U - w/z)B i 2 > b > (-t 4 ) + 9(y - z - t 3 )B ( $ (-t 3 ) 



X s 



0(ta - u,/z)B«(-i 3 )} + 9(4™ - 1) {l?«>(t§) + B^(tS) + + ^HS)} + 9(1 - 4™)9(y - z) x 

B<$(w/z) - B 2 j b \-w/z) ) - [9(1 - 4w)9(4w - y 2 ) - 9(y 2 - Aw)] Q(z-y) x 



,0) ' 



.CO 



7T\/ — 



jy2 



2s 



CO 



' B^iw/z) - B^{-w/z) 



(B30) 



1 

-Rex 



(i) 



2B< J c ) (f 1 ) - 9(i 2 - w/z)B${w/z) - S(w/z - t 2 )B^{t 2 ) - ®{h - z + y)Bl J c >(t 1 ) 



Mi 



-9(y - z - t^B^i-h) + e(z + y- h)B 2 J c >(-z - y) + 9(ti - z - y)B^J{-t 1 ) + <d{w/z - t 2 )B^{-t 2 ) 
-2Bg\-t 2 ) - 29(y 2 - Aw) [B? c \t 3 ) + B^(t 4 )] - 2S(Aw - y 2 ) [b«(^) + B«(t|)] + 9(1 - Aw) x 

x {9(i 4 - w/z)B$(w/z) + Q(w/z - U)B [ ^{U) + 0(t 3 - z + y)B^(t 3 ) + 9(y - z - t 3 )B^(-t 3 ) 

+[0(i 4 - z - y) + 0(t 3 - z - y)]B ( £{-z - y) - Q(w/z - U)0(U - z ~ y)B^(-U) 
- 0(t 3 - z - y)B«(-t 3 )} + 9(4 W - 1) {B^{t%) + B^{tt) + B&\-t c 3 ) + B&\-tl)} 

+9(1 - 4w)9(y - z) (s< j) ™ +s^' ) 7T-24 7 ' ) [b^(w/z) - B^(-z - y)]) - 9(1 - 4w)9(z - y) x 



;(.;) 



(J), 



M n V z2 -y 2 



- 2s 



0) 



B^Kz)-BW(-z-y) 



(B31) 



1 



.0) 



2S«(ix) - B«(ix) - B«(-tx) - 2B { i(z - y) + 2B$(-t 2 ) - 2B«(t 4 ) 



- Re x d 



+^(t 3 )+^(-t 3 )+2S^(f 3 )+ S 



B«(u,/z)-B^(z-y) 



(B32) 
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In these expressions we have used along with (IB27[) 
and (|B28|) the following arguments 



\/y 2 +Aw =F ; 



-y t vV 



*i,a 

^3,4 

t% A = ±[i T zV^T^T; 

l 3,4 



4w 



(B33) 
(B34) 
(B35) 
(B36) 



Note that i§ 4 and i§ 4 are complex- valued, which as- 
sumes the analytic continuation of elliptic functions 
hinted in Appendix B of Ref. 



Transforming Re% 



terms of and Bif ' , we accumulate the residue terms 
which are different for each response function. They are 
fully defined by the following sets of constants 



into the representation in 



S « = 
e? = 

,(3,4) 
'2 

S (5) _ 



,( 5 ) - 1 



(1) 



,(1) 



4 2) = 4 2) = o, 



„(3,4) 



;(3.4) 



,(5) 



(B37) 
(B38) 
(B39) 
(B40) 



We recall once again that for j ' = 3 and j = 5 
Eqs. (|B29]) - (|B32|) must be complemented by x (j) ([BlT]> . 



1 S. A. Wolf ef aZ., Science 294, 1488 (2001). 

2 I. Zutic, J. Fabian, and S. Das Sarma, Rev. Mod. Phys. 
76, 323 (2004). 

3 Y. A. Bychkov and E. I. Rashba, J. Phys. C 17, 6039 
(1984). 

4 J. Nitta, T. Akazaki, H. Takayanagi, and T. Enoki, Phys. 
Rev. Lett. 78, 1335 (1997). 

5 T. Koga, J. Nitta, T. Akazaki, and H. Takayanagi, Phys. 
Rev. Lett. 89, 046801 (2002). 

6 E. G. Mishchenko and B. I. Halperin, Phys. Rev. B 68, 
045317 (2003). 

7 E. G. Mishchenko, A. V. Shytov, and B. I. Halperin, Phys. 
Rev. Lett. 93, 226602 (2004). 

8 A. Burkov, A. S. Nunez, and A. H. MacDonald, Phys. Rev. 
B 70, 155308 (2004). 

9 A. G. Mal'shukov, L. Y. Wang, C. S. Chu, and K. A. Chao, 
Phys. Rev. Lett. 95, 146601 (2005). 

10 R. Raimondi, C. Gorini, P. Schwab, and M. Dzierzawa, 
Phys. Rev. B 74, 035340 (2006). 

11 M. Pletyukhov, Phys. Rev. B 75, 155335 (2007). 

12 A. Shnirman and I. Martin, Europhys. Lett. 78, 27001 
(2007). 

13 L.I. Magarill, A.V. Chaplik, and M.V. Entin, JETP 92, 
153 (2001). 

14 E. Rashba, Phys. Rev. B 70, 161201(R) (2004). 

S. I. Erlingsson, J. Schliemann, and D. Loss, Phys. Rev. B 
71, 035319 (2005). 

16 J. A. Maytorena, C. Lopez-Bastidas, and F. Mireles, Phys. 
Rev. B 74, 235313 (2006). 

17 A. Brataas, A. G. Mal'shukov, and K. A. Chao, Phys. Rev. 
B 56, R1684 (1997). 



18 C. A. Ullrich and M. E. Flatte, Phys. Rev. B 66, 205305 

(2002) . 

19 W. Xu, Appl. Phys. Lett. 82, 724 (2003). 

20 X. F. Wang, Phys. Rev. B 72, 085317 (2005). 

21 G. Gumbs, Phys. Rev. B 72, 165351 (2005). 

22 M. S. Kushwaha and S. E. Ulloa, Phys. Rev. B 73, 205306 
(2006). 

23 M. Pletyukhov and V. Gritsev, Phys. Rev. B 74, 045307 
(2006). 

24 F. Stern, Phys. Rev. Lett. 18, 546 (1967). 

25 H. Imamura, P. Bruno, and Y. Utsumi, Phys. Rev. B 69, 
121303(R) (2004). 

26 M. A. Ruderman and C. Kittel, Phys. Rev. 96, 99 (1954). 

27 T. Kasuya, Prog. Theor. Phys. 16, 45 (1956). 

28 K. Yosida, Phys. Rev. 106, 893 (1957). 

29 G. D. Mahan, Many-Particle Physics (Plenum Press, New 
York, 1990). 

30 I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Se- 
ries, and Products (Academic, New York, 1994) . 

31 V. M. Edelstein, Solid State Commun. 73, 233 (1990). 

32 E. I. Rashba and Al. L. Efros, Appl. Phys. Lett. 83, 5295 

(2003) . 

33 A. Shekhter, M. Khodas, and A. M. Finkel'stein, Phys. 
Rev. B 71, 165329 (2005). 

34 D. C. Hamilton and A. L. McWhorter, in Light Scattering 
Spectra of Solids, edited by G. B. Wright (Springer, New 
York, 1969). 

35 Guang-Hong Chen and M. E. Raikh, Phys. Rev. B 60, 
4826 (1999). 



